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With a radio continuum galaxy survey by the Square Kilometre Array (SKA), a photometric 
galaxy survey by Euclid and their combination, we forecast future constraints on primordial non- 
Gaussianity. We focus on the potential impact of local-type higher-order nonlinear parameters on 
the parameter estimation and particularly the confirmation of the inflationary consistency inequality. 
Nonstandard inflationary models, such as multifield models, introduce the scale-dependent stochastic 
clustering of galaxies on large scales, which is a unique probe of mechanism for generating primordial 
density fluctuations. Our Fisher matrix analysis indicates that a deep and wide survey provided by 
SKA is more advantageous to constrain tnl, while Euclid has a strong constraining power for /nl 
due to the redshift information, suggesting that the joint analysis between them is quite essential 
to break the degeneracy between the nonlinear parameters. The combination of the full SKA and 
Euclid will achieve the precision level needed to confirm the consistency inequality even for /nl ~ 1-5 
and Tnl ~ 17, though it is still hard for a single survey to confirm it when /nl fS 2.7. 

PACS numbers: 


I. INTRODUCTION 

Observing the large-scale distribution of matter pro¬ 
vides us the rich information about not only the late-time 
evolution but also the primordial nature of the Universe. 
Among various cosmological parameters characterizing 
the primordial Universe, we are particularly interested in 
possible departures from a purely Gaussian distribution 
of primordial density fluctuations called primordial non- 
Gaussianity (PNG), which is one of the most powerful 
tests of inflation and more generally a key to understand¬ 
ing the extreme high-energy physics. 

One of the major theoretical discoveries is that all in¬ 
flationary models predict the presence of the consistency 
relation between the parameters characterizing PNG. For 
the simplest scenarios, if there is the nonvanishing local- 
form bispectrum, the trispectrum must necessarily exist 
with Tnl = ((6/5)/nl)^, where /nl and tnl are the so- 
called local-type nonlinear parameters characterizing the 
amplitude of the primordial bispectrum and trispectrum. 
Even in a general situation, one can show that there is 
a universal relation, tnl > ((6/5)/nl)^, which is often 
refereed to as the Suyama-Yamaguchi inequality 
(see also [1, |3|)- The primordial trispectrum is also usu¬ 
ally characterized by another nonlinear parameter, ^nl, 
which corresponds to the strength of the intrinsic cubic 
nonlinearities of primordial fluctuations. Hence a detec¬ 
tion of the higher-order PNG and the confirmation of 
the inequality would indicate the presence of more com¬ 
plicated dynamics, e.g., multifield inflationary models, in 
the primordial Universe and should be, thus, the target 
in future experiments (see also S). 


‘Email: yamauchi”at”jindai.jp 


The current limits on these parameters have been 
obtained from cosmic microwave background (CMB) 
anisotropies: = 0.8 ± 5.0 and = 

(—9.0 ± 7.7) X 10^ at Icr statistical significance p and 
^ 2800 (95% G.L.) 0. However, CMB measure¬ 
ments are already close to being cosmic-variance lim¬ 
ited. A complementary way to access PNG is to mea¬ 
sure its impact on a large-scale structure. For the stan¬ 
dard Gaussian initial conditions, the halo bias is often 
assumed to be linear, deterministic, and scale indepen¬ 
dent. It was found that the non-Gaussianity in primor¬ 
dial fluctuations effectively introduces a scale-dependent 
clustering of galaxies on large scales iQ. In this pa¬ 
per, we focus on scale-dependent stoehastic halo bias. If 
"T^L > ((6/5)/nl)^, the halo density contrast is not 
100% correlated to the matter linear density field 5 even 
in the absence of the shot noise ©[niiii. Namely, the 
simple bias relation dh = bh^ should be modified due to 
the stochasticity. The halo-halo power spectrum between 
b- and 6'-mass bins, cannot be expressed as the 

products of the linear bias defined by = Pj^l/Pmm 
and the linear matter density power spectrum Pnim- For¬ 
mally this means P^’^ ^ > b^^b^ ^Pmm- This provides 
the unique opportunity of large-scale scale-dependent 
stochastic bias as a probe of the primordial Universe as¬ 
sociated with complicated dynamics. 

A promising way to explore PNG in future large-scale 
structure surveys is to reduce the cosmic-variance noise 
with the so-called multitracer technique [l3,13 : which 
allows us to measure the scale-dependent bias very accu¬ 
rately. Since using the multitracer technique we can con¬ 
strain the PNG in the halo bias without suffering from 
the cosmic-variance noise, the clustering analysis in this 
case is expected to be limited only due to the contri¬ 
butions from the shot noise. In order to reduce the re¬ 
maining noise contribution, the two-dimensional map of 
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the large-scale structure projected along the line of sight 
is considered, while the redshift information of galaxies 
would be lost by the projection. Although the multi¬ 
tracer technique is truly effective for three dimensional 
statistics, in this paper we will focus on two-dimensional 
one as a simple extension of our previous work M- Even 
when we use the projected density contrasts as tracers, 
the multitracer technique is expected to be still effec¬ 
tive. Indeed Ref. [l^ shows that the constraining power 
on /nl drastically improves even when splitting samples 
into two tracers and increases with the number of trac¬ 
ers. In this paper we extend the multitracer technique to 
include the effects of the higher-order PNG and discuss 
the required survey level needed to test the consistency 
relation. As future representative surveys, we consider 
Euclid- ^ in optical and infrared bands, and the Square 
Kilometre Array (SKA)- ^ in radio wavelength. Both 
telescopes will perform ultimate galaxy surveys with dif¬ 
ferent characteristics. The Euclid photometric galaxy 
survey (15,000 deg^) will reach z fv 2.5 and provides red- 
shift information, which is highly advantageous to con¬ 
strain /nl- On the other hand, the SKA continuum 
galaxy survey will cover a wider area of sky (30,000 deg^) 
and a significant redshift depth {z < 5), while the red- 
shift information is not available. As we will show below, 
the wide and deep survey provided by SKA can constrain 
tnl effectively. 


II. ANGULAR POWER SPECTRUM 

In what follows, we compute the angular power spec¬ 
trum for the halo density contrast, extending our previ¬ 
ous work [T^ to include the effects of tnl and ^nl (see 


also [laiii)- We first split whole samples into mass- 
divided subsamples for each redshift bin to apply the 
multitracer technique. It can be shown that in the pres¬ 
ence of /nl and tnl the halo-halo power spectrum be¬ 
tween the b th and b' th mass bins is expressed as [isini 
(see also El) 
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where = j5^p/D+{z) is the growth factor, 
and M{k) = 2fc^r(fc)/3nni,o-ffo ^ith T(k) being the 
matter transfer function normalized to unity at large 
scales El- Throughout the paper, we use the expres¬ 
sion = 25c{b^Q — 1), where 6c is the critical density 
for spherical collapse and we will take the <5c = 1.46 to fit 
the numerical simulation [^ . We employ a fit to simu¬ 
lation for the mass function dn/dM and the linear bias 
factor b^Q given in [2l|. If there is a nonvanishing ^nl, 
we should take into account the additional correction to 
the halo bias as EH) 1^ - 

and we adopt the fitting function for to simula¬ 
tion [^. When defining /nl, tnl, and (?nl, we have 
evaluated primordial fluctuations at present, though at 
the decoupling for the CMB convention, suggesting that 
the observed nonlinear parameters have the relation [ 2 I : 
/nl 1.3/^^®, Tnl ~ 1.32t^^®, and ^nl ~ 1.3ggf®. 
With Eq. Ill), the halo-halo angular power specta be¬ 
tween b- and 5'-mass bins in the i th redshift bin is given 
by 
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where d^V/dzdD, = x^/H, x is the coinoving distance, 
and = fo°° denotes the av¬ 

erage density. We have introduced to represent the 
selection function. 

To take advantage of the multitracer technique, we 
need to estimate the halo mass of each galaxy, which 
has to be inferred from available observables. How¬ 
ever, since estimates of the halo mass involve large un¬ 
certainties, a number of nuisance parameters should be 
included to model systematic errors. We assume the 
mass-observable relation including uncertainties in the 
mass inference from available data. In our treatment, 
the probability of assigning the estimated mass Mggt 
to the true mass M is assumed to be given by log¬ 
normal distribution with the variance crinM(AI, z) and 


the bias lnMbias(Af, z) [ 2 ^. With these, the selection 
function can be expressed as 5'i(b)(M, z) = r(t,)0(z — 
Zi)0(zi+i-z)i [erfc (a;(A/({,); M))-erfc (a;(M(b+i); A/))], 
where a;(Mest; M) = (InMest-lnM-lnMbiasVv^crinM, 
r(f,) is the gray-body factor to denote the ratio between 
halos and what we really observe for each mass bin, be¬ 
cause we may not observe all galaxies associated with 
the underlying dark matter halos. Hereafter we will in¬ 
troduce 14 nuisance parameters in the variance and the 
bias to quantify the impact of possible residual system¬ 
atic errors on the parameter estimation [l^ . . 

In order to understand the constraints on the non- 
Gaussian parameters, /nLi^nl and 5 nl, it is useful to 
note that their effects on the halo-halo power spec- 
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trum have different redshift-, scale- and halo-mass de¬ 
pendences. The correction to the bias roughly scales as 
Pfz/k'^ for /nl, for tnl, and j5gz/k'^ for ^nl- 

Because the correction to the bias from /nl is identi¬ 
cal to the one from ^nl = (/3///3g)/NL, /nl and ^nl 
are degenerate for a single tracer case. However, in the 
case with multiple tracers, the different halo-mass depen¬ 
dence of/?/ and fig would break their degeneracy [l^. On 
the other hand, because the bias correction due to tnl 
has larger dependences on the scale and redshift, its de¬ 
tectability would be enhanced for a survey with larger 
angular scale and higher redshift coverage. Thus, we ex¬ 
pect that, with the multitracer technique, wide and deep 
surveys will be powerful to probe higher-order PNG. 


III. RESULTS 

To study the required survey level needed to test 
the consistency relation, we proceed to the Fisher 
analysis. The Fisher matrix is defined by Fap = 

where the indices /, J 
run over the redshift and mass bins (/,6,6'), 0“ are 
free parameters to be determined by observations, and 
Ci{(.) = is the observed power 

spectrum including the shot noise contamination. The 
marginalized expected Icr error on 0“ is estimated to be 
ct(0“) = \/ {F~^)aa- We adopt the covariant matrix for 
multiple tracers whch are observed in different sky areas 
with some overlap given in [l^ . 

Before showing the expected constraints, we need to 
specify the survey parameters. We consider the SKA 
radio continuum survey, the Euclid photometric survey, 
and their combination as future representative surveys 
with significant high precisions. SKA covers 30,000deg^ 
out to 2 = 5, though there is only one redshift bin since 
the redshift information is not available. Then we sim¬ 
ply drop the redshift dependence in and InMSKA, 

To infer the halo masses, we consider five radio galaxy 
types such as star-forming galaxies, radio quite quasars, 
radio-loud AGN (FRI and FRII), and starbursts. Fol¬ 
lowing [s^, we will assign these galaxy types the fol¬ 
lowing halo mass: {A/sfG: ?Wrqq, Mfri, Msb, AfpRii} = 
{1,30,10^,5 X 10^,10^} in the unit of In 

order to have a plausible distribution for the halo 
mass associated with each population, we introduce 
the five separating masses M(i) = 0.9 x IO^^/i^^Mq, 

M(2) = y^MsfG^RQQ, Af(3) = \/-^^RQQ-^FRI, Af(4) = 

VA/friA/sb, and Af( 5 ) = V^'^SbAAfrii- With these, 
the five mass bins can be defined through AL(j) < 
M < AL(j+i) (i = 1,2,3,4) and M > To 

match the expected number density distribution of galax¬ 
ies [ 23 , [13: adopt the gray-body factor = 

{0.013,0.03,0.1,1,1} and ^ {0.2,0.4,1,1,1}, re¬ 

spectively. As for Euclid, the covered area 15,000deg^ 
and the redshift range 0.2 < z < 2.7 are considered. The 



FIG. 1: Forecast Icr marginalized contours in (/nliLnl), 
(gNLiRNL), and (/nl,(?nl) planes with the vanishing fiducial 
values of the nonlinear parameters. To see the dependence on 
the minimal multipole we also plot the results for SKA with 
= 3 in the dashed line. 


redshift information is provided though photometric red- 
shifts. Galaxy samples are then split into five redshift 
bins with same interval Az = 0.5. Since photometric 
surveys provide various galaxy properties such as lumi¬ 
nosity, color, and stellar mass, which can be used to infer 
the halo mass, we can further split the galaxy samples ac¬ 
cording the estimated halo mass. We consider five mass 
bins such that each mass bin of the same redshift bin has 
the same number density, presumably because the tight¬ 
est constraint is expected to be obtained when the shot 
noises for the mass bins become comparable. For the 
flux cut, we adopt the minimum observed mass for each 
redshift bin, AAest.min/(IO^^/ i'^ALq) = 0.7,1,2, 5,10, and 
set = 1 instead. For the combination of these 

surveys, the area of the overlap region is assumed to 
be 9,000deg^ and we neglect the contributions from the 
derivative of the cross-correlations for simplicity. In total 
we include 22 = 8(SKA)-|-I4(Euclid) nuisance parame¬ 
ters to model the systematic errors as well as three non¬ 
linear parameters. We choose = ^.3, 

and zero for the other parameters as fiducial values. To 
calculate the Fisher matrix, we use 2 < ^ < 400 for 
SKA and 3 < £ < 400 for Euclid. We hereafter fo¬ 
cus on constraints on /nl,t-nl, and ^nl, marginalizing 
over the other parameters. Our fiducial model is a stan¬ 
dard ACDM cosmological model with the parameters: 

= 0.318, Hb,o = 0.0495, Oar = 0.6817, w = -1, 
h = 0.67, rig = 0.9619, /cq = 0.05Mpc“^, and erg = 0.835. 

In order to see the impact of the higher-order PNG 
on the parameter estimation, we first plot the Icr ex¬ 
pected marginalized contours in Fig. [U in the case with 
vanishing nonlinear parameters. Although the resultant 
constraints on /nl are slightly weaker than the previous 
results [I1I13 where tnl and ^nl were neglected, the 
constraints from both SKA and Euclid are still signifi- 
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FIG. 2: Parameter space to confirm the consistency relation, 
/NL/cr(/NL) > 1 and T>jL/o'('rNL) > 1. For comparison, the 
inconsistent region and the constraints from Planck are also 
shown in gray color. 


stochastic value, there is little stochasticity and the tight¬ 
est constraint on /nl is obtained, as expected. On the 
other hand, with increasing tnl, the constraining power 
on /nl decreases, mainly because the correction from tnl 
to the halo bias dominates that from /nl- Hence, the rel¬ 
atively smaller tnl is needed to detect /nl- We find from 
the figure that even for SKAl there is a wedge-shaped re¬ 
gion where we can confirm the consistency inequality at 
more than Icr level. As is anticipated, the confirmation of 
the inequality becomes harder as decreasing /nl- When 
/nl is small, say < 2.7, the confirmation at the > Icr level 
is rather challenging for a single survey, even with the 
multitracer technique. However, combining SKA2 and 
Euclid can break the degeneracy between /nl and tnl, 
and drastically improve the constraints to reach the con¬ 
firmation of the consistency relation even for /nl ~ 1.5 
and Tnl ~ 17. We roughly estimate the boundary of the 
viable region as /nl ^ 1 and 10 < tnl ^ 1-5(/nl)^- 


IV. SUMMARY 


cant, cr(/NL) = 0{1). Especially, the redshift informa¬ 
tion obtained from the photometric survey significantly 
improves the constraint on /nl 0- As for tnl, SKA2 
can reach cr(rNL) = 23, which is an improvement by a 
factor of 100 compared with the Planck constraint. In 
contrast, the constraint from Euclid, ^(tnl) = 62, is rel¬ 
atively weaker than one from SKAl, ^(tnl) = 43. This is 
understood as follows. The bias correction from tnl has 
stronger dependence on the wavelength and redshift than 
those from /nl and 5 nl as argued above. Then wider sky 
coverage and deeper redshift are expected to be more ad¬ 
vantageous to constrain tnl - Although the redshift infor¬ 
mation is not available for SKA, the contributions from 
high-z samples are significant to the effective bias. To 
see the effect of sky coverage, we also plot the expected 
contours for SKA with -^min = 3 in the dashed line. We 
find that ^rain = 3 reduces the power of the SKA signif¬ 
icantly and becomes comparable or weaker than Euclid, 
suggesting that the full-sky observation of the SKA is es¬ 
sential when we try to constrain the higher-order PNG. 
Similarly, the expected constraints on ^nl are given as 
cr(gNL) = 1-8 X 10^ (Euclid), 2.1 X 10"^ (SKAl), and 
7.4 X 10^ (SKA2), which are several tens of times sev¬ 
erer than that obtained from Planck. These different 
behaviors suggest that SKA and Euclid turned out to be 
complementary probes of PNG and the joint analysis be¬ 
tween SKA and Euclid are quite effective to confirm the 
PNG consistency inequality. 

We then study the dependence of our forecast on the 
choice of the fiducial values, concentrating on /nl and 
Tnl- For this purpose, we consider stochastic bias with 

values of Tnl which satisfies Tnl > ((6/5)/nl)^- In Fig. [21 

we show the region where both /nl and tnl are de¬ 
tected at ItT significance, that is, /nl/o’(/nl) > 1 and 
tnl/o-(tnl) > 1 in (/nLi^nl) plane for SKA1(2), Euclid, 
and their combinations. When tnl is close to the non- 


To summarize, we have discussed the potential impact 
of the higher-order nonlinear parameters tnl )5nl, and 
the required survey level needed to test the PNG consis¬ 
tency relation tnl > ((6/5)/nl)^ for future galaxy sur¬ 
veys such as SKA radio continuum and Euclid photomet¬ 
ric surveys. With the multitracer technique, the Fisher 
matrix analysis has revealed that deep survey and large 
sky coverage provided by SKA are advantageous to con¬ 
strain Tnl, while the Euclid photometric survey has the 
strong constraining power for /nl thanks to redshift in¬ 
formation. The information from both Euclid and SKA is 
quite essential to break the degeneracy between the non¬ 
linear parameters. Indeed the combination of SKA2 and 
Euclid can detect the consistency inequality in the wide 
parameter region at more than Icr level, more specifically 
even for /nl ~ 1-5 and tnl ~ 17, though for a single sur¬ 
vey it is still hard to confirm when /nl ^ 2.7. 

Our analysis also implies that the large-angle obser¬ 
vations is quite essential to provide the constraints on 
the higher-order PNG. We should address the spectro¬ 
scopic surveys conducted by both the SKA and Euclid. 
In this paper we conservatively assumed no redshift in¬ 
formation for the SKA and relatively large redshift bin 
for the Euclid. The redshift information for individual 
galaxies obtained from the spectropic surveys may allow 
the tomographic analysis to enhance the Fisher matrix 
due to the cross-correlations between the different red¬ 
shift bins, providing the improvement of the constraints. 
We hope to come back to these issues in the near future. 

When the constraints on nonlinear parameters, in par¬ 
ticular /nl, are close to 0(1), the general relativistic 
correction in the observed power spectrum might not be 
negligible [3ll-[^. We have simply neglected this effect 
expecting that we can subtract the effect from the ob¬ 
served power spectrum. 
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